Chemical and mechanical instability in warm and dense nuclear matter 
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We investigate the possible thermodynamic instability in a warm and dense nuclear medium 
(r < 50 MeV and po < ps < 3po) where a phase transition from nucleonic matter to resonance- 
dominated A matter can take place. The analysis is performed by requiring the global conservation 
of baryon and electric charge numbers in the framework of a relativistic equation of state. Similarly 
to the liquid-gas phase transition, we show that the nucleon-A matter phase transition is character- 
ized by both mechanical instability (fluctuations on the baryon density) and by chemical-diffusive 
instability (fluctuations on the charge concentration) in asymmetric nuclear matter. We then per- 
form an investigation and a comparative study on the different nature of such instabilities and phase 
transitions. 

PACS numbers: 21.65.-f, 25.75.-q, 64.10.-|-h 



I. INTRODUCTION 

One of the most interesting aspects of the experiments 
on heavy-ion collisions is a detailed study of the thermo- 
dynamical properties of strongly interacting nuclear mat- 
ter away from the nuclear ground state. In this direction, 
many efforts have been focused on searching for possible 
phase transitions in such collisions. At low temperatures 
(T < 10 MeV) and subnuclear densities, a liquid-gas type 
of phase transition was first predicted theoretically 
and later observed experimentally in a nuclear multifrag- 
mentation phenomenon at intermediate-energy nuclear 
reactions i5j. 

Because nuclei are made of neutrons and protons, the 
nuclear liquid-gas phase transition is in a binary system 
where one has to deal with two independent proton and 
neutron chemical potentials for baryon number and elec- 
tric charge conservation. Taking into account of this im- 
portant property, a very detailed study of Miiller and 
Serot ^ focused on the main thermodynamic properties 
of asymmetric nuclear matter in the framework of a rel- 
ativistic mean field model. 

A relevant aspect of a system with two conserved 
charges (baryon and isospin numbers) is that the phase 
transition is of second order from the viewpoint of Ehren- 
fests definition. At variance with the so-called Maxwell 
construction for one conserved charge, the pressure is not 
constant in the mixed phase and therefore the incom- 
pressibility does not vanish 0| • Such feature plays a 
crucial role in the structure and in the possible hadron- 
quark phase transition in compact star objects [1, Q. 
Moreover, for a binary system with two phases, the bin- 
odal coexistence surface is two dimensional and the in- 
stabilities in the mixed liquid-gas phase arise from fluctu- 
ations in the proton concentration (chemical instability) 
and in the baryon density (mechanical instability) [lol - 

[a. 

Although the equation of state (EOS) at densities be- 
low the saturation nuclear matter is relatively well known 
due to the large amount of experimental nuclear data 



available, at larger densities there are many uncertain- 
ties; the strong repulsion at short distances of nuclear 
force makes, in fact, the compression of nuclear matter 
quite difhcult. However, in relativistic heavy-ion colli- 
sions the baryon density can reach values of a few times 
the saturation nuclear density po and/or high tempera- 
tures. The future CBM (compressed baryonic matter) 
experiment of the FAIR project at GSI Darmstadt will 
make it possible to create compressed baryonic matter 
with a high net baryon density [l^. In this direction 
very interesting results have been obtained at low energy 
at the CERN Super Proton CoUider (SPS) and at low- 
energy scan at the BNL Relativistic Heavy Ion Collider 

(RHIC) dMl]. 

In regime of finite values of density and temperature, 
a state of high density resonance matter may be formed 
and the A(1232)-isobar degrees of freedom are expected 
to play a central role in relativistic heavy ion collisions 
and in the physics of compact stars |17H21j . Transport 
model calculations and experimental results indicate that 
an excited state of baryonic matter is dominated by the 
A-resonance at the energy from the BNL Alternating 
Gradient Synchrotron (AGS) to RHIC [H-il. More- 
over, in symmetric nuclear matter and in the framework 
of a non-linear Walecka model, it has been predicted that 
a phase transition from nucleonic matter to A-excited 
nuclear matter can take place and the occurrence of this 
transition sensibly depends on the A-meson coupling con- 
stants [1^, [131 • Due to the presence of only one conserved 
"charge" (baryon number) considered in these previous 
investigations, the region of the phase transition develops 
when the incompressibility becomes negative and there- 
fore only mechanical instabilities are present. 

The information coming from experiments with heavy 
ions in intermediate- and high-energy collisions is that 
the EOS depends on the energy beam but also sensibly on 
the electric charge fraction Z/A of the colliding nuclei, es- 
pecially at not too high temperatures 28, 29]. Moreover, 
the study of nuclear matter with arbitrary electric charge 
fraction results to be important in radioactive beam ex- 
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periments and in the physics of compact stars. 

In this article, we study the hadronic EOS at finite 
temperature and density by means of a relativistic mean- 
field model with the inclusion A-isobars and by requiring 
the Gibbs conditions on the global conservation of baryon 
number and net electric charge. In this context, let us 
observe that, for the range of temperatures and baryon 
densities considered in this investigation (T < 50 and 
PB < 3/90 ), the contribution of strange hadron particles 
can be neglected in a good approximation due to their 
very low concentration. In fact, unlike compact stars in 
a /3-stability regime, since weak decays cannot take place 
during the short lifetime of a high density system, the 
only possibility of producing strangeness is through asso- 
ciated production but, in the scenario we are discussing, 
this process has been shown to be very inefficient [3^, iSlJ 
and, therefore, the study of the possible phase transition 
can be limited to two conserved charges. 

The main goal of this paper is to show that, for an 
asymmetric warm and dense nuclear medium, the pos- 
sible A-matter phase transition is characterized by me- 
chanical and chemical-diffusive instabilities. Similarly to 
the liquid-gas phase transition, chemical instabilities play 
a crucial role in the characterization of the phase transi- 
tion and can imply a very different electric charge fraction 
Z / A in the coexisting phases during the phase transition. 

The paper is organized as follows. In Sec. II, we 
present the relativistic hadronic equation of state. In 
Sec. Ill, we review the most important thermodynamic 
proprieties of the phase transitions in a binary system, 
highlighting the relevant features of mechanical and dif- 
fusive instabilities. The main results are presented in 
Sec. IV, which is divided into two subsections: in A, we 
review the most important results obtained in the liquid- 
gas phase transition and, in B, we study the nucleon-A 
matter phase transition. Finally, in Sec. V, we summa- 
rize our conclusions. 



II. HADRONIC EQUATION OF STATE 

The relativistic mean-field model (RMF) is widely suc- 
cessful used for describing the properties of finite nuclei 
as well as hot and dense nuclear matter (32| - [36| . 

In the RMF model the Lagrangian density for nucleons 
can be written as 

-Cn = Vat 7m ^ i^'^N - 9a N (y) - guiN 7m "^^ 

-9pN 7m t- P'] ^JV + \{d^,(jd^(T - mla^) - U{a) 

+ \ mi wpw^ + i c {gl^ oj^ajf")^ ^'^ ' 

-\f,,f>^'^ ~\g,M^\ (1) 

where Mj^ = 939 MeV is the nucleon vacuum mass and 
t is the isospin operator which acts on the nucleon. The 
field strength tensors for the vector mesons are given by 



the usual expressions F^i, = d^uj^, — dyW^, G^i, = d^pi, — 
dvP^n and U{a) is the nonlinear potential of a meson 

U{a) = ^a{gaNaf + -^b^g^Nof , (2) 

usually introduced to achieve a reasonable compression 
modulus for equilibrium normal nuclear matter [ssj . In 
the following, the meson-nucleon coupling constants and 
the other parameters (a, 6, c) of the EOS will be fixed to 
the parameters set marked as TMl of Ref. [35| . 

In a regime of finite values of temperature and den- 
sity, a state of high-density resonance matter may be 
formed and the A(1232)-isobar de gree s of freedom are 
expected to play a central role [itI llSl. [2^ [23 | . In par- 
ticular, the formation of resonances matter contributes 
essentially to baryon stopping, hadronic fiow effects and 
enhanced strangeness [37|. 

It is well known that thus far there is no relativistic 
quantum theory for the A as a spin 3/2 field without 
any inconsistency when imposing other fields such as the 
ones with electromagnetic interaction [s^. Moreover, fol- 
lowing the Rarita-Schwinger formalism, the spin 3/2 par- 
ticle, described by means of a vector spinor state, has an 
off-shell spin 1/2 sector. To incorporate A-isobars in the 
framework of effective hadron field theories, a formalism 
was developed to treat A analogously to the nucleon, tak- 
ing only the on-shell As into account and the mass of the 
As are substituted by the effective one in the relativistic 
mean field approximation [39l . |40[ . The Lagrangian den- 
sity concerning the A-isobars can be then expressed as 
27, 40, 41] 

^A=i>Ai' [*7m^^ ^ (^^A - .go-AO-) - gu,Alfi^^'']i>A J (3) 

where V'a Rarita-Schwinger spinor for the A- 

isobars (A++, A"*", A'^, A~). Due to the uncertainty 
on the meson-A coupling constants, we limit ourselves to 
consider only the coupling with the a and w meson fields, 
more of which are explored in the literature (see Sec. IV 
for details) 113, 111142]. 

In the RMF approach baryons are considered as Dirac 
quasiparticles moving in classical meson fields and the 
field operators are replaced by their expectation values. 
As a consequence, the field equations in a mean field 
approximation are 





- M*j^ - guiNl^t^ - gpNl'^hp)i)N 


= 0, (4) 




-Ml-g^A7M^A = 0, 


(5) 


mla + 


«5^ArO-^ + bglj^a^ = ^ g„,pf , 

i 


(6) 


mluj + 


i 


(7) 


mlp^ 




(8) 



i 



where a = {a), uj = (a;°) and p = {p^}, are the non- 
vanishing expectation values of meson fields, the index i 
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runs over the considered baryon particles, and the effec- 
tive mass of the zth baryon is defined as 



(9) 



The pf and pf are the baryon density and the baryon 
scalar density, respectively. They are given by 



1: 



d^fc 



[n^(fc) - ni(k)] , 



(27r)3 ' 



(10) 

(11) 



where 7^ is the degeneracy spin factor (7JV = 2 and 
7A = 4) and ni{k) and ni{k) are the fermion particle 
and antiparticle distribution functions, given by 



ni{k) = 
ni{k) = 



1 



exp(i?*(fc)-/i*)/T+l' 
1 

exp(i?*(fc) + M*)/r + l- 



(12) 
(13) 



The baryon effective energy is defined as Ei*{k) — 

fc^ + Mi*^ and the effective chemical potentials p,* are 
given in terms of the meson fields as follows: 



(14) 



where pi are the thermodynamical chemical potentials. 
Pi = de/dpi. 

Because we are going to describe a finite temperature 
and density asymmetric nuclear matter, we have to re- 
quire the conservation of two " charges" : baryon number 
[B) and electric charge (C) (as already remarked, we 
neglect the contribution of strange hadrons, because a 
tiny amount of strangeness can be produced in the range 
of temperature and density explored in this study). As 
a consequence, the system is described by two indepen- 
dent chemical potentials: pB and pc, the baryon and the 
electric charge chemical potential, respectively. There- 
fore, the chemical potential of particle of index i can be 
written as 



Pi = bi PB +CiPC : 



(15) 



where bi and Ci are, respectively, the baryon and the 
electric charge quantum numbers of the ith hadron. 

The thermodynamical quantities can be obtained from 
the baryon grand potential 57^ in the standard way. More 
explicitly, the baryon pressure Pb = —^Ib/V and the 



energy density can be written as 



Pb = 



IE 



^22 
1 22 



1 2 2 



U{a) 



[ni{k) + ni{k)] 



^22 
2 " 



NUJ) 



(16) 



d^fc 
(2^ 



E*{k) [n,(fc) + n,(fc)] 



1 



Uia) + -m^uj + - c {gujN w) 



(17) 



Let us observe that the contribution of the lightest non- 
strange mesons (pious) may not be negligible in regime 
of temperature and density achieved during the possible 
A-matter phase transition. Following Ref. [i^, from a 
phenomenological point of view, we can take into account 
the pion degrees of freedom by adding their one-body 
contribution to the thermodynamical potential, that is, 
the contribution of an ideal Bose gas with an effective 
pion chemical potential p*, depending self-consistently 
from the meson fields. The value of /i* is obtained from 
the "bare" one Ptt, given from Eq. p^ . and subse- 
quently expressed in terms of the corresponding effective 
baryon chemical potentials, respecting the strong interac- 
tion. More explicitly, fromEq. P3|) . — pc = Pp~Pn 
and the corresponding effective chemical potential can be 
written as 

= /4 - f^*n^ f^p- ~ 9pN P, (18) 

where the last equivalence follows from Eq. As a 

consequence, = and /Lt*o = p-„o = 0. 

This assumption can be seen somehow in analogy with 
the hadron resonance gas within the excluded-volume 
approximation. There the hadronic system is still re- 
garded as an ideal gas but in the volume reduced by 
the volume occupied by constituents (usually assumed 
as a phenomenological model parameter), here we have 
a (quasifree) pion gas but with an effective chemical po- 
tential that contains the self-consistent interaction of the 
meson fields. 

Of course, this naive phenomenological approach can- 
not incorporate the very complex TriVA interaction at fi- 
nite temperature and baryon density and a more realistic 
chiral symmetric model should be implemented. On the 
other hand, as we will see in Sec. IV-B, such an effective 
nuclear EOS has the noticeable advantage of simplify- 
ing the not trivial numerical analysis involved in seeking 
of thermodynamic instabilities and in the construction of 
the mixed phase. Due to this fact, it would be prudent to 
see the results of this preliminary study in a perspective 
of academic interest. 

Finally, the total pressure and energy density are given 
by the baryon (B) and pion (M) contribution: e — es + 
EM and P = Pb + Pm- 
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III. PHASE TRANSITIONS AND STABILITY 
CONDITIONS 

As already stated, we are dealing with the study of a 
multi-component system at finite temperature and den- 
sity with two conserved charges: baryon number and elec- 
tric charge. For such a system, the Helmholtz free energy 
density F can be written as 

F{T, pB,Pc) = -P{T, PB, Pc) + PbPb + P-cPc , (19) 



with 



fJ-B = 



( dF 



\9pi 



P-c = 



T,pc 



dF 



(20) 



In a system with N different particles, the particle chem- 
ical potentials are expressed as the linear combination of 
the two independent chemical potentials p b and pc and, 
as a consequence, X]il=i MiPi = PbPb + PcPc- Therefore, 
the number of particles may change during a process and, 
at variance of density and temperature, different particle 
degrees of freedom may be relevant in the description of 
the system (for example, at low temperature and den- 
sity, we have protons and neutrons only, while at higher 
temperature and density other kind of particles, such as 
A-isobars, can appear). What it is actually relevant for 
the thermodynamical description under consideration are 
only the two conserved charges and not the number of 
different particles constituent the system. 

In general, a system can exist in a number of different 
phases, each of which exhibit quite different macroscopic 
behavior. The single phase that is realized for a given 
set of independent variables is the one with the lowest 
free energy. In a system with two conserved charges, it 
is possible to have iVmax = 4 phase coexistence regions 
in thermodynamical equilibrium [4^ l45j , even if we have 
found no evidence for the existence of more than two 
phases in the regime investigated in this paper. There- 
fore, assuming the presence of two phases (denoted as 
/ and //, respectively), the system is stable against the 
separation in two phases if the free energy of a single 
phase is lower than the free energy in all two phases con- 
figuration. The phase coexistence is given by the Gibbs 
conditions 



f^B 



P\T,pb,pc)^P"{T,pb,pc)- 



(21) 
(22) 



higher density (//) phase and x is the volume fraction of 
the phase // in the mixed phase (0 < x < 1)- 

An important feature of this conditions is that, unlike 
the case of a single conserved charge, the pressure in the 
mixed phase is not constant and, although the total pb 
and Pc are fixed, baryon and charge densities can differ 
in the two phases, according to Eq.s ([23|) and ((24|). 

For such a system in thermal equilibrium, the possi- 
ble phase transition can be characterized by mechanical 
(fiuctuations in the baryon density) and chemical insta- 
bilities (fiuctuations in the electric charge density). As 
usual the condition of the mechanical stability implies 



PB 



/ dP 
\dpB 



> 0. 



(25) 



T, PC 



By introducing the notation p^ j = (dpi/ dpj)T,p (with 
i, j = B,C), the chemical stability can be expressed with 
the following conditions 45] 



PB,B > , pc,c > , 



Pb,b Pb,c 
PC.B pc,c 



>0. (26) 



In addition to the above conditions, for a process at con- 
stant P and T, it is always satisfied that 



Pb Pb,b + PC Pc,B = , 
Pb Pb,c + PC Pc,c = . 



(27) 
(28) 



Whenever the above stability conditions are not re- 
spected, the system becomes unstable and the phase 
transition take place. The coexistence line of a system 
with one conserved charge becomes in this case a two di- 
mensional surface in (T, P, y) space, enclosing the region 
where mechanical and diffusive instabilities occur. 



IV. RESULTS AND DISCUSSION 

As mentioned in the Introduction, our main goal is 
to study the instability regions related to the formation 
of the A-isobars at finite temperature and baryon den- 
sity. Because such instabilities may have several analo- 
gies with the liquid-gas phase transition, it is instructive, 
first, to briefly review the main properties of this nu- 
clear phase transition in the framework of our EOS and 
to test the numerical procedure that will be applied at 
higher temperatures and densities. 



Therefore, at a given baryon density pB and at a given 
net electric charge density pc — y Pb (with y — Z/A), the 
chemical potentials p b are pc are univocally determined 
by the following equations 

PB = (1-x)Pb{T.,Pb-,Pc)+XP^{T,pb,Pc)-, (23) 
PC = (1 - x) PciT,PB,Pc) + XPc{T,PB,Pc) , (24) 



A. Liquid-gas phase transition 

In a regime of low temperature and baryon density, 
relevant in the liquid-gas phase transition, only proton 
and neutron degrees of freedom take place. In this simple 
case, for example, Eq. ((28)) can be written as 



where p^ 



and Pc^^^ are, respectively, the baryon and 
electric charge densities in the low density (/) and in the 



dpp 
dy 



= 0. 



(29) 



T,P 



5 



P [MeV/fm^] 




-0.2^ 



FIG. 1: (Color online) Pressure as a function of baryon den- 
sity for various values of the proton fraction. The continuous 
(dashed) lines correspond to the solution obtained with (with- 
out) the Gibbs construction. 



[MeV] 




FIG. 2: (Color online) Proton and neutron chemical poten- 
tial as a function of the proton fraction y for various isobars 
(P=0.25, 0.20, 0.15, 0.10, 0.075 MeV/fm^) (lines a to e) at 
r=10 MeV. 



where y — pp/ ps- Because we are working with a proton 
fraction < y < 0.5, the chemical stability conditions 
(1^5)) are therefore satisfied if 

( ^2/ )t,p ^ ^ (a?/)y p'^^ ^^'^^ 

[due the validity of Eq. ([2^ . the first above condition 
implies the second one and vice versa] . 

As aheady observed, in presence of two conserved 
charges the liquid-gas phase transition can be charac- 
terized by mechanical and chemical instabilities Q. In 
order to better put this feature in focus, we report in 
Fig. [1] the pressure as a function of baryon density for 
various values of the electric charge fraction y at fixed 
temperature T = 10 MeV. The continuous lines corre- 
spond to the solution obtained with the Gibbs construc- 
tion, whereas the dashed lines are without correction. 
For a proton fraction y > 0.2 a mechanical instability is 
present, whereas for y < 0.2 the system becomes unstable 
only under chemical-diffusive instability. 

The presence of chemical unstable regions are much 
more evident in Fig. [5J where we show the proton and 
neutron chemical potentials for various isobars at con- 
stant temperature, as a function of the proton asymme- 
try. Below P — 0.25 MeV/fm^, the system becomes 
unstable because of the presence of regions of negative 
(positive) slope for pp {p,„). 

In order to study the phase coexistence of the system, 
in Fig. [3l we show the binodal section as a function of 
the proton asymmetry y at T = 10 MeV. Following the 
same notation of Ref. \ff\ , the binodal surface is divided 
into two branches by a critical point (CP) and a point 
of equal equilibrium (EQ) at = 0.5, where protons and 
neutrons have the same concentration. The left branch of 
the diagram represents the initial phase configuration of 
the system at lower density (gas phase, I) and the second 
branch, at higher density, corresponds to the final phase 
configuration (liquid phase, II). 

The binodal surface encloses the area where the system 
undergoes to the phase transition. The mixed-phase re- 



gion extends up to small values of the proton asymmetry, 
whereas the mechanical instability region ends around 
y^O.2. 

During the isothermal compression, the system evolves 
through configuration at constant y and meets the first 
branch in a point A. At this point the system becomes 
unstable and an infinitesimal phase in B appears at the 
same temperature and pressure of A. In this context 
let us remember that, although the proton asymmetry is 
globally conserved, this is not true for the single phase. 
In particular, for an asymmetric nuclear system it is ener- 
getically favorable to separate it into a liquid phase (less 
asymmetric) and a gas phase (more asymmetric) rather 
than into two phases with equal proton fraction. 

If point A has a value of yA greater than the corre- 
sponding values ycp of the CP (as in the upper panel of 
Fig. [3]), the system ends the phase transition in the liquid 
phase (in the point C). On the other hand, as already 
observed in Ref. @ , if the system has been prepared in 
a very asymmetric configuration with ?/^/ < ycp (lower 
panel of Fig. it undergoes to a retrograde phase tran- 
sition. A second liquid phase in B is formed but, after 
reaching a point of maximum volume fraction Xmax < 1, 
the system returns to its initial gas phase at point C . 
Note that this kind of phase transition is possible only 
for a multi-component system and in this case, is purely 
diffusive. 

To better characterize the evolution of the two phases, 
in Fig. m the volume fraction x of the second phase 
during the phase transition is shown. By increasing the 
asymmetry parameter of the system under consideration 
(at lower values of y), the maximum density achieved 
during the mixed phase decreases, until the system un- 
dergoes a retrograde phase transition (for y < 0.15). 
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FIG. 3: (Color online) Binodal section at T = 10 MeV, with 
in evidence the critical point (CP), the point of maximum 
asymmetry (MA) and the point of equal equilibrium (EQ). In 
the upper and lower panels are reported the evolution of the 
mixed phase for two different system configurations (see the 
text for details). 
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FIG. 4: (Color online) Evolution of the volume fraction x of 
the second phase as a function of the baryon density for a 
system with different values of j/ at T = 10 MeV. 



B. A-matter phase transition 

By increasing the temperature and the baryon density 
during the high energy heavy ion collisions (T « 50 MeV 
and po < Pb < Spo), a multi-particle system with A- 
isobar and pion degrees of freedom may take place. 




Pb/ Po 



FIG. 5: The energy per baryon versus baryon density at zero 
temperature and y — 0.5 with (a) no A; (b) x^a = 1.10; (c) 
x^A = 1.15; (d) x^A = 1.20; (e) x^a = 1.27; (f) XaA = 1.30; 
(g) XcrA = 1.33. 



To better understand the relevance of A-isobar and the 
dependence of the EOS on the meson-A coupling con- 
stants {XaA = 9aA/gcrN, X^iA = Swa/SljTv), in Fig. O WC 

report the energy per baryon as a function of the baryon 
density at zero temperature and y — 0.5 for different val- 
ues of XaA and x^jA = 1- Let note that by increasing 
the value of x^Ai a second minimum on the energy per 
baryon appears. 

FoUovifing Ref. [4l| . in setting x^a and x^^a, we have 
to require that (i) the second minimum of the energy per 
baryon lies above the saturation energy of normal nu- 
clear matter, i.e., in the mixed A-nucleon matter only a 
metastable state can occur; (ii) there are no A-isobars 
present at the saturation density; and (iii) the scalar 
field is more (equal) attractive and the vector potential is 
less (equal) repulsive for As than for nucleons, in accor- 
dance with QCD finite-density calculations 42]. In this 
context, it is proper to remember that QCD sum-rule 
predictions for the scalar self-energy are sensitive to the 
unknown density dependence of four-quark condensates 
and, due to this, there is no certainly reliable information 
about the coupling constant of the A isobars with scalar 
mesons. 

Of course, the choice of couplings that satisfy the above 
conditions is not unique but exists a finite range of pos- 
sible values (represented as a triangle region in the plane 
XcA-Xuja) which depends on the particular EOS under 
consideration [4l[. Without loss of generality, in the fol- 
lowing we can limit our investigation to move only in 
a side of such a triangle region by fixing XuiA = 1 and 
varying x^a h'om unity to a maximum value compati- 
ble with the aforementioned conditions. As can be ob- 
served in Fig. [51 for the TMl parameter set, such a max- 
imum value corresponds to = 1.33, while the value 
x^^^ = 1.27 corresponds to the appearance of the second 
minimum on the energy per baryon with the formation 
of a metastable state. Analog behaviors can be obtained 
with other EOS parameters set (see, for example, Ref. 
'^I'l for more details). 

In Fig. [SI we show in symmetric nuclear matter the 
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FIG. 6: (Color online) The relative nucleons (solid lines) and 
A (dashed lines) densities as a function of the baryon density 
for different values of temperature with (a) a^o-A ~ 1-2; (b) 
XaA = 1-27; (c) XaA — 1-33. 



in the range 1 < x^a < 1-1, instabilities are restricted to 
very low values of temperature and electric charge frac- 
tion, but for Xa-A > 1-1, such instabilities start to be 
much more relevant and extend to higher values of T 
and y. 



P [MeV/fm^] 




Pfl/PO 



relative nucleon (solid lines) and the A-isobar (dashed 
lines) density fraction {Yi = pi/pB) versus the baryon 
density at T = and T = 50 McV, for different values 
of XaA- We observe that A-matter becomes dominant 
with respect to the nucleon concentration at high baryon 
density and such effect is significantly anticipated by in- 
creasing the temperature [iq. 

In analogy with the liquid-gas case, we are going to 
investigate the existence of a possible phase transition in 
the nuclear medium by studying the presence of instabil- 
ities (mechanical and/or chemical) in the system. 

As already observed, during a phase transition with 
two conserved charges, the electric charge fraction y = 
Pel Pb is not locally conserved in the single phase but 
only globally conserved. Therefore, during the compres- 
sion of the system, the appearance of particles with nega- 
tive electric charge (such as A~) could, in principle, shift 
the diffusive instability region to negative values of ?/, 
even if the system is prepared with a positive y. Such a 
feature has no counterpart in the liquid-gas phase tran- 
sition and, as we will see, it turns out to be very rele- 
vant in order to properly determine the instability region 
through the binodal phase diagram. 

Taking into account that Eq. becomes in this case 



dy 



T,P 



C 



dy 



= 0. 



(31) 



T,P 



the chemical stability condition is satisfied if 



dy 



{(^) <0,if,>0, 



> or 



T,P 



( dp.1 



. \ dy 



> , if ?/ < . 



T.P 



(32) 

It is relevant to observe that for the value x^a = 1, we 
do not find any mechanical or diffusive instability. Con- 
trariwise, by increasing the x^a coupling ratio, mechan- 
ical and chemical instabilities take place. In particular. 
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FIG. 7: (Color online) Pressure as a function of baryon 
density for T = 30 MeV (upper panel) and T = 50 MeV 
(lower panel) with x^a = 1.3. Labels a to f correspond to 
y = 0.5, 0.4, 0.3, 0.2, 0.1, 0, respectively. 



To better clarify this aspect, we report in Fig. [71 the 
pressure as a function of the baryon density at T = 30 
MeV (upper panel) and T = 50 MeV (lower panel), for 
different values of y and Xo-a = 1-3. In this context it 
is interesting to observe that at T = 50 MeV and below 
y = 0.3, the system becomes mechanically stable, but, 
in a similar manner to the liquid-gas case, is chemically 
unstable. This important feature can be better observed 
in the Fig. |S1 where we report the baryon and electric 
charge chemical potential isobars as a function of y, at 
fixed temperature T = 50 MeV and Xcta = 1-3. 

From the analysis of the above chemical potential iso- 
bars, we are able to construct the binodal surface relative 
to the nucleon-A matter phase transition. In Fig. |9l we 
show the binodal section at T = 50 MeV and XaA = 1-3. 

The right branch (at lower density) corresponds to the 
initial phase (I), where the dominant component of the 
system is given by nucleons. The left branch (II) is re- 
lated to the final phase at higher densities, where the 
system is composed primarily by A-isobar degrees of free- 
dom (A-dominant phase). In the presence of A isobars. 
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FIG. 8: (Color online) Baryon (upper panel) and electric 
charge (lower panel) chemical potential isobars as a function 
of 2/ at T = 50 MeV and x^a = 1-3. The curves labeled 
a through g have pressure P—9, 7, 6, 5, 4, 3, 2 MeV/fm'^, 
respectively. 



the phase coexistence region results differ substantially 
from what was obtained in the liquid-gas case, in partic- 
ular it extends up to regions of negative electric charge 
fraction and the mixed phase region ends in a point of 
maximum asymmetry with y — —1 (corresponding to a 
system with almost all A~ particles, with the contribu- 
tion from antiparticles and pions almost negligible in this 
regime) . 

Repeating the reasoning made for the liquid-gas phase 
transition, we analyze the phase evolution of the system 
during the isothermal compression from an arbitrary ini- 
tial point A, indicated in Fig. [S] In this point, the system 
becomes unstable and starts to be energetically favorable 
the separation into two phases; therefore, an infinitesimal 
A-dominant phase appears in _B, at the same tempera- 
ture and pressure. Let us observe that, although in B the 
electric charge fraction is substantially negative, the rel- 
ative A~ abundance must be weighed on the low volume 
fraction occupied by phase II near point i? (x ~ 0). Dur- 
ing the phase transition (0 < x < l)j each phase evolves 
towards a configuration with increasing y, in contrast to 
the liquid-gas case, where each phase evolves through a 
configuration with a decreasing value of y (with the ex- 
ception of the gas phase after the maximum asymmetry 
point). 
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FIG. 9: (Color online) Binodal section at T = 50 MeV and 

XcrA = 1.3. 
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FIG. 10: (Color online) Volume fraction of the A-matter 
phase as a function of the baryon density for a system with 
different values of y. 



In order to better understand the evolution of the two 
phases in the mixed phase, we report in Fig. [10] the 
volume fraction x of the A-matter phase as a function 
of the baryon density. Unlike in the liquid-gas case, by 
decreasing the electric charge fraction y of the system 
under consideration, the mixed phase involves a greater 
region of baryon density and extends below the nuclear 
saturation density. 

In the previous example we considered a fixed value of 
temperature. The maximum temperature at which the 
system becomes mechanically stable depends on the par- 
ticular value of the electric charge fraction. For example, 
at y = 0.3, it is about T,„ax = 49.5 MeV and at y = 0.5, it 
is about Tmax — 50.6 MeV. Furthermore, when y = 0.5, 
the end of the mechanical instability region, obviously, 
corresponds to the end of the mixed-phase region. This 
is no longer true in the presence of two conserved charges. 
Due to the presence of a diffusive instability region, the 
mixed phase can extend to slightly higher temperatures 
with respect to the maximum temperature achieved in 
the symmetric case. Although this feature involves small 
differences in temperature, is interesting from a concep- 
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FIG. 11: (Color online) Pressure as a function of baryon den- 
sity at different values of y, from y = 0.5 (label a) to y = 
(label f). The continuous (dashed) lines correspond to the 
solution obtained with (without) the Gibbs construction. 



tual point of view to investigate this aspect in more de- 
tail. Toward this purpose, in Fig. [11] we show the pres- 
sure as a function of the baryon density and the Gibbs 
construction (continuous lines) for various values of the 
electric charge fraction, x^a — 1-3 and T = 51 MeV 
(dashed lines are without Gibbs construction). In this 
case the system is always mechanically stable, while it is 
unstable for the presence of the chemical-diffusive insta- 
bility up to y = 0.35. 

At lower temperatures, the mixed-phase region be- 
comes more relevant at higher values of y. This feature 
can be seen in Fig. 1121 where the binodal section (upper 
panel) and the isothermal pressure as a function of the 
baryon density (lower panel) is reported at T = 40 MeV 
and icrA = 1-3. The Gibbs construction corresponds to 
the curve from A to C; the isothermal curves in B and 
D (with ub 7^ vd) are also reported. In this case, we 
assume that the system is initially prepared in the low- 
density (nucleonic) phase with y = 0.3, corresponding 
to point A. During the compression each phase evolves 
following the corresponding curve up to points C and 
D, where the system leaves the instability region in the 
A-matter phase. 

In Fig. [131 we show the Gibbs construction (continuous 
lines) to the EOS aX y — 0.3, XcrA — 1-3 and for different 
temperatures. By decreasing the temperature, the insta- 
bility region extends over a wide range of baryon density. 
In particular, below T = 40 MeV, the phase transition 
starts slightly below the nuclear saturation density. 

As already observed, the mixed phase structure results 
strongly affected not only by the temperature, but also 
by the particular choice of the XaA coupling. In fact, by 
decreasing the a- A coupling constant, the mixed phase 
region shifts to lower temperatures. To better clarify this 
aspect, we study the phase transition for x^-a = 1-22 and 
T = 20 MeV (at T = 50 MeV, the system resuhs to be 
mechanically and chemically stable). 

In Fig. 1141 upper panel, we report the pressure as 
a function of the baryon density for different values of 
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FIG. 12: (Color online) Binodal section at r= 40 MeV and 
XcA = 1.3, with the point of equal equilibrium in evidence. 
(Lower panel) The corresponding isothermal curves, with the 
Gibbs construction (curve from the point A to C) at y = 0.3 
and the isotherms of points B and D shown. 
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FIG. 13: (Color online) Isotherms at constant y = 0.3 and 
Xo-A = 1.3, for various values of temperatures. The solid 
(dashed) lines represent the EOS obtained with (without) 
Gibbs construction. 



y. The continuous lines correspond to the Gibbs con- 
struction in the region of instability of the EOS. For this 
choice of parameters, the binodal section (lower panel) 
differs substantially with respect to the previous cases 
and two separate regions of instability are present. The 
first one extends at lower pressure and it is present only 
for small value of y, where both mechanical and diffu- 
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FIG. 14: (Color online) (Upper panel) Pressure as a function 
of the baryon density for different values of y from y — 0.5 
(label a) to y = (label f), with the instability regions (Gibbs 
construction in the continuous lines) in evidence. (Lower 
panel) The binodal diagram, with the point of equal equi- 
librium and the two critical points in evidence. In both insta- 
bility sectors a region of retrograde phase transition is present. 



sive instabilities are present. Let us observe that in this 
lower region, for y > ycp, the system goes to a retro- 
grade phase transition likewise to the liquid-gas phase 
transition. The upper region of instability extends at 
greater pressures and higher values of y, where mechan- 
ical and diffusive instabilities are both present. Also in 
this second region, on the left of the CP, a retrograde 
phase transition can occur. However, in this particular 
case the system is already in a A-dominant phase and, 
at the end of the mixed phase, in which A-isobars are 
partially converted into nucleons, it quickly returns to 
the A-matter phase. 

Finally, in Fig. I15[ we report the phase diagram with 
in evidence the coexistence regions of the liquid-gas and 
the nucleon-A matter phase transition for y = 0.3 and 0.5 
{xaA — 1.3). The two coexistence regions are well sepa- 
rated and the features of the two phase transitions differ 
significantly. In fact, for the liquid- gas transition, asym- 
metric nuclear matter implies a reduction of the second 
critical density and of the critical temperature Tc- Con- 
trariwise, for the A-dominant phase transition, we have 
a slight increase of the critical temperature and a signif- 
icant reduction of the first critical density. In particular 
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FIG. 15: (Color online) Phase diagram of the liquid-gas and 
the nucleon-A matter phase transition for y = 0.3 (dashed 
curves) and y — 0.5 (continuous curves). The lines labeled 
with I and II, delimitate the first and second critical densities 
of the coexistence regions, respectively. 



at moderate temperatures (T « 30 ~ 40 MeV), the sys- 
tem begins the mixed phase at a baryon density of the 
order of pQ. This behavior could be phenomenologically 
relevant in order to identify such a phase transition in 
heavy ion collision experiments. 

V. CONCLUSIONS 

The main goal of this work is to show the possible ex- 
istence of chemical and mechanical instability at finite 
temperature and dense nuclear matter. We have stud- 
ied the relativistic nuclear EOS with the inclusion of A- 
isobars and by requiring global conservation of baryon 
and electric charge numbers. Similarly to the liquid-gas 
phase transition in a warm and low density nuclear mat- 
ter, a nucleon-A matter phase transition also can oc- 
cur at higher temperatures and densities (T < 50 MeV, 
p « 1 ^ 3po)- We have shown that for asymmetric nu- 
clear matter both mechanical and chemical instabilities 
take place. The latter plays a crucial role in the char- 
acterization of the phase transition and can also imply 
very low values of the electric charge fraction y during 
the mixed phase region. 

The nucleon-A matter phase transition depends signif- 
icantly on the value of the ct-A coupling constant and we 
have seen that only for a limited range of the possible 
(physical) x^a couplings the presence of thermodynamic 
instabilities may become relevant from a phenomenolog- 
ical point of view. 

Whether metastable A-excited nuclear matter exists is 
still a controversial issue because little is actually known 
about the A-coupling constants with the scalar and vec- 
tor mesons, even if QCD finite-density sum rule results 
predict a larger net attraction for a A-isobar than for a 
nucleon in the nuclear medium (4^ . Although we have 
seen that instabilities are already present for Xo-a > 1, 
they become phenomenologically more relevant at greater 
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values of Xo-a , involving larger region of mixed phase and 
greater values of the electric charge fraction. 

The analysis of the instability regions with different 
A-coupling constants turns out to be not trivial from 
the numerical point of view, especially at lower values of 
XctA where a complex structure of the mixed phase can 
be formed. For example, we have shown that, in the case 
of Xo-A = 1.22, two separate regions of instability are 
present. Moreover, the case XaA = 1-3 has been studied 
for different values of temperature and we have seen that 
in asymmetric nuclear matter the mixed phase transition 
involves a large range of baryon densities. 

Similarly to the liquid-gas phase transition, the nu- 
cleonic and the A-matter phase have a different electric 
charge fraction in the mixed phase. The electric charge 
fraction in the nucleonic phase reflects a system with 
higher values of y than the A-matter phase. In the liquid- 
gas phase transition, the process of producing a larger 
neutron excess in the gas phase is referred to as isospin 
fractionation |Tl|, |l3] ■ A similar effects can occur in 
the nucleon-A matter phase transition due essentially to 
a A~ excess in the A-matter phase with lower values of 
y. As already observed, due to the uncertainty on the 
meson-A coupling constants, we have not considered in 
this investigation the coupling of the A with the isovector 
p-meson field, because this is much less explored in the 
literature. We have verified that the presence of a p — A 
coupling could further increase the isospin asymmetry in 
the mixed phase and lower the critical temperature of the 
nucleon-A matter phase transition. 

In this context it is proper to observe that Coulomb 
interaction and finite size effects, not considered in this 
study, can significantly alter the structure of the phase 
transition. Moreover, as already observed, it is proper to 



remember that our effective EOS cannot incorporate the 
complex ttA^ A dynamics and it would be very interesting 
to investigate the presence of thermodynamic instabili- 
ties in the framework of a more realistic chiral hadronic 
EOS. Taking also into account the large uncertainty on 
the possible values of A-meson field couplings, it would 
be prudent to note the pedagogical character of this pre- 
liminary study. 

Many effects discussed in this paper may be more ev- 
ident at low values of y, obtainable, in principle, with 
radioactive ion beam facilities. On the other hand, it is 
rather unlikely, at least in the near future, that neutron 
rich nuclei can be accelerated to energies larger than a 
few GeV per nucleon. However, some precursor signals 
of the considered instabilities could be observed even in 
collisions of stable nuclei at intermediate energies. For 
example, in Ref. [2^, the simulation of the reaction 
238u_^238u (average y = 0.39), at 1 A GeV and semi- 
central impact parameter 5 = 7 fm, shows that rather 
exotic nuclear matter can be formed in a transient time 
of the order of 10 fm/c, with a baryon density up to 
3 po, < 50 H- 60 MeV and y « 0.35 0.40. Such condi- 
tions would agree fully with the results for the nucleon-A 
mixed phase region (see Fig. [T5|) . 

A possible signature of the nucleon-A matter phase 
transition could be found via observables particularly 
sensitive to the expected different isospin content of the 
two phases. For example, at the AGS energies, the A- 
resonance was predicted to be the dominant source for 
pious of small transverse momenta [13] ■ In this case, an 
increase of the negative pious 7r~ of small trasverse mo- 
menta at a greater asymmetry of the beam could be a 
good indicator of a A isospin fractionation effect. 
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